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1 Introduction and overview 

Living cells are soft bodies of a characteristic form, but endowed with a 
capacity for a steady turnover of their structures. Both of these material 
properties, i.e. recovery of the shape after an external stress has been im- 
posed and dynamic structural reorganization, are essential for many cellular 
phenomena. Examples are active intracellular transport, cell growth and 
division, and directed movement of cells. The structural element responsi- 
ble for the extraordinary mechanical and dynamical properties of eukaryotic 
cells is a three-dimensional assembly of protein fibers, the cy to skeleton. A 
major contribution to its mechanical properties is due to actin filaments 
and proteins that crosslink them. Numerous experiments in vivo ( [Eichinger 



3t al., 1996| ) and in vitro (Janmey, 1991; Janmey, 1995| ) have shown that 



the mechanical properties of cells are largely determined by the cytoskeletal 
network. 

The cytoskeletal polymers (actin filaments, microtubules and intermedi- 
ate filaments) which build up this network are at the relevant length-scales 
(a few microns at most) all semiflexible polymers. In contrast to fiexible 
polymers the persistence length of these polymers is of the same order of 
magnitude as their total contour length. Hence, the statistical mechanics 
of such macromolecules can not be understood from the conformational en- 
tropy alone but depends crucially on the bending stiffness of the filament. 



* To be published in "Dynamical networks in physics and biology" (Proceedings of a Les- 
Houches workshop), edited by G. Forgacs and D. Bey sens, EDP Sciences Springer Verlag 
(1998). 
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The nontrivial elastic response and distribution function of the end-to-end 
distance of a single semiflexible polymer illustrate instructively this inter- 
play of energetic and entropic contributions. Fig. || shows the extension (in 
units of fLp'/ksT) of a semiflexible polymer of fixed length L grafted at one 
end when a weak force / is applied at the other end for different persistence 
lengths Ip ( |Kroy and Frey, 1996 ). This is one of the few exact results for 
the wormlike chain model. Note that the largest response is obtained for a 
polymer with a persistence length of the order of its contour length. More 
flexible molecules contract due to the larger entropy of crumpled conforma- 
tions, stiffer molecules straighten out to keep their bending energy low. In 
the flexible case the response is isotropic and proportional to l/ksT, i.e., 
the Hookian force coefficient is proportional to the temperature and we re- 
cover ordinary rubber elasticity. On the other hand, when the persistence 
length is longer than the contour length, the response becomes increasingly 
anisotropic. Transverse forces give rise to mechanical bending of the fila- 
ments and the transverse spring coefficient in the stiff limit is proportional 
to the bending modulus k. The effective longitudinal spring coefficient turns 
out to be proportional to k^/T, indicating the breakdown of linear response 
at low temperatures T ^ (or ip — > oo). This is a consequence of the well 
known Euler buckling instability illustrated in Fig. |l] (top). The anisotropy 
of the elastic response considerably complicates the construction of network 
models. 

Energetic effects also have a strong influence on the probability distribu- 
tion function G{R; L) for the end-to-end distance i? of a semiflexible poly- 
mer of length L. When the polymer is not much longer than its persistence 
length, G{R] L) deviates strongly from the Gaussian shape found for flexible 
polymers (see Fig. |2|) with the main weight being shifted towards full exten- 
sion with increasing stiffness ( [Wilhelm and Frey, 1996| ) . This is also reflected 
by the nonlinear force-extension relation, whose overall shape is determined 
by the Euler instability. Measurements of the radial distribution function 
as well as of the response of single polymers to applied forces can be used 
to measure the bending stiffness of biopolymers. 

A useful experimental technique for investigating the short time dynam- 
ics of semiflexible polymers is dynamic light scattering ( |Farge and Maggs 



1993| ; Harnau et al., 1996| ; Kroy and Frey, 1997 ). We have calculated the 



dynamic structure factor for solutions of semiflexible polymers in the weakly 
bending rod limit. The result allows characteristic parameters of the poly- 
mers, such as their persistence length ip, their lateral diameter a and also 
the mesh size of the network to be determined from light scattering mea- 
surements. The most important result of this work is that the decay of the 
dynamic structure factor shows two regimes, an initial simple exponential 
decay determined by the hydrodynamics of the solvent and a stretched expo- 
nential decay due to structural relaxation driven by the bending energy. For 
actin and intermediate filaments, where both regimes can be resolved, the 
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Figure 1: Top: The elastic response of a stiff rod is extremely anisotropic due to 
the Euler instability. Bottom: Response of a filament clamped at one end with a 
fixed initial orientation to a small external force at the other end. The normalized 
extension (inverse force coefficient) is plotted versus the persistence length in units 
of the total length of the polymer. The curves are parameterized by the relative 
angle 9o of the force f to the initial orientation. 



method is well suited to investigate also more complex questions concerning 
for example the interactions with actin binding proteins which among others 
can induce cross-linking or bundling. 

In order to describe the material properties of the cytoskeleton, one has 
to understand how semiflexible polymers built up statistical networks and 
how the macroscopic stresses and strains are mediated to the single filaments 
in such a "rigid polymer network". This is a quite rapidly growing field 



in polymer physics with many questions still open (Aharoni and Edwards, 
T99|) 



If the polymers are not crosslinked, the response will depend on how fast 
one pulls. Roughly speaking, the solution will either show elastic behavior 
and obey Hook's law with a linear relation a = G'j between stress a and 
strain 7 or rather show viscous behavior and obey Newton's law a = rjj, 
where the stress is proportional to the strain rate 7. The task of a theoretical 
description is to find out how the material parameters like the shear modulus 
and the viscosity r] depend on the elastic and dynamic properties of 
the polymers and the architecture of the network. The mechanical and 
dynamical properties of the various crosslinking proteins are also expected 
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to influence the viscoelasticity of the network (Wachsstock et al., 1994). 

Whereas it is known that on short time scales the effect of mutual steric 
hindrance ( "entanglement" ) in conventional polymeric materials is very sim- 
ilar to the effect of permanent chemical crosslinks, this is an open problem 
for semiflexible polymer networks. Since forces between neighboring poly- 
mers can only be transmitted transverse to the polymer axis and there is 
no restoring force for sliding of one filament past another, a steric contact is 
acting completely different from a chemical crosslink on a microscopic scale. 

But even the most simple of these problems, the static elasticity of a 
network of permanently crosslinked semiflexible polymers, is quite complex. 
Because of the strongly anisotropic behavior of the single elements, the pre- 
dicted macroscopic properties of the network vary greatly with the explicit or 
implicit assumptions made about network geometry and stress propagation 
by recent theoretical treatments of both entangled solutions and crosslinked 
networks ([MacKintosh et al., 1995; [sambert and Maggs, 199(3; Satcher and 



Dewey, 1996| ; [Kroy and Frey, 19961 ). We hope to represent a key aspect of 
the geometrical structure of both cellular and artificial stiff polymer net- 
works by looking at disordered networks. Specifically, we use a crosslinked 
network of sticks randomly placed in a plane as a toy model for studying 
the origin of macroscopic elasticity in a stiff polymer network. Although 
quantitative predictions about the behavior of existing (three-dimensional) 
networks of semiflexible polymers are not attempted at this stage, this model 
is expected to reflect the salient features of the full problem and to pro- 
mote its understanding by allowing the detailed discussion of questions like 
"What modes of deformation contribute most to the network elasticity?", 
"How many filaments do actually carry stress, how many remain mostly 
unstressed?", "What kind of effective description of the complicated mi- 
croscopic network geometry should be used?". This approach connects the 
theory of cytoskeletal elasticity to the very active fields of transport in ran- 
dom media and elastic percolation. In section 3^ we will describe the model 
in more detail and discuss the question of the dominant deformation mode. 

Previous fluorescence microscopic observations ( |Kas et al., 1994| ) suggest 
that the tube picture is a useful concept for understanding the viscoelastic 
properties of entangled semiflexible polymer networks. Starting from this 
model, where the matrix surrounding a single test polymer is represented 
by a tube-like cage, we have developed a phenomenological description that 
seems to be able to account for some of the observed elastic and dynamic 
properties. The predictions for the entanglement transition, the concentra- 
tion dependence of the plateau modulus, and the terminal relaxation time 
are in good quantitative agreement with experimental data. 



2 Single-chain properties 
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2.1 The wormlike chain model 

The theoretical understanding of the mechanical properties of a single semi- 
flexible macromolecule in isolation is already a nontrivial statistical mechan- 
ics problem with quite a number of recent developments 50 years after it was 
first formulated ( [Kratky and Porod, 194^ ). The model usually adopted for 
a theoretical description of semiflexible chains like actin filaments is the 
wormlike chain model. The filament is represented by an inextensible space 
curve r(s) of total length L parameterized in terms of the arc length s. 
The statistical properties of the wormlike chain are determined by a free 
energy functional which measures the total elastic energy of a particular 
conformation. 

n= I ds^f?^]'; |t| = l, (1) 



ds) 

where t(s) = dr{s)/ds is the tangent vector. The energy functional is 
quadratic in the local curvature with k being the bending stiffness of the 
chain. The inextensibility of the chain is expressed by the local constraint, 
|t(s)| = 1, which leads to non-Gaussian path integrals. Since there would 
be high energetic costs for a chain to fold back onto itself one may safely 
neglect self-avoidance effects for sufficiently stiff chains. 

Despite the mathematical difficulty of the model some quantities can be 
calculated exactly. Among these is the tangent-tangent correlation function 
which decays exponentially, (t(s) • t(s')) = exp [— (s — s')/ip], with the per- 
sistence length ip = K/ksT (in three dimensional space). Another example 
is the mean-square end-to-end distance 

7^2 := (R^) = 2ll{e-^'^^ - 1 + L/ip) 

for L/£p ^ (rigid rod) 
for L/ip oo (random coil), 

which reduces to the appropriate limits of a rigid rod and a random coil (with 
Kuhn length 2ip) as the ratio of L to ip tends to zero or infinity, respectively. 
The calculation of higher moments quickly gets very troublesome ( Herman^ 
and Ullman, 1952| ). 




2.2 Linear force-extension relation 

Another useful property of the model, which can be computed exactly is the 
linear force extension relation ( [Kroy and Frey, 1996 ). Adding a term — f • R 



to the Hamiltonian in Eq. [T|, the extension of a wormlike chain with a weak 
force applied between its ends is computed as 

6R = {{R)f- (R)o) f // = fiTZ^ - (|R|)')/A:,r = f/k (3) 
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with k the force coefficient. The problem is the calculation of the moment 
of the the end-to-end distance (|R|). This has not been achieved so far, 
but an expansion in the stiff limit is possible. For example, we can write 
(|R|) = 7^(y^1 + — TZ'^)fR?) and expand the square root to obtain to 
leading order for the force coefficient 

4A;bT- 



This is exact in the stiff limit, L/lp <C 1, and qualitatively correct over 
the whole range of stiffness. However, for the special boundary condition 
of a grafted chain we can even get an exact result for arbitrary stiffness. 
Consider a chain with one end clamped at a fixed orientation and let us 
apply a force at the other end. Then the linear response of the chain may 
be characterized in terms of an effective Hookian spring constant which 
depends on the orientation Oq of the force with respect to the tangent vector 
at the clamped end. In the appropriate generalization of the second part 
of Eq. ^ the force coefficient k is replaced by an angle dependent function, 
k~^ — > Lp'Tof^^lp/ L) /ksT . Noting that the conformational statistics of the 



worm like chain is equivalent to the diffus ion on the unit sphere ( Saito et al. 



1967| ) the function T can be calculated ( Kroy and Frey, 1996|) (see Fig 



In the flexible limit, where the chain becomes an isotropic random coil, 
all curves coincide and reproduce entropy elasticity, Teo{x) ~ x. But for 
stiff chains the force-extension relation depends strongly on the value of the 
angle 9q between the force and the grafted end. Transverse forces give rise 
to ordinary mechanical bending characterized by the bending modulus k 
~ 1/2^)) whereas longitudinal deformations are resisted by a larger 
force coefficient k? /T {!Fq ~ l/x^). The breakdown of linear response in the 
limit T ^ (or £p 00) is a consequence of the well known Euler buckling 
instability, as we already mentioned in the introduction. The linear response 
for longitudinal forces is due to the presence of thermal undulations, which 
tilt parts of the polymer contour with respect to the force direction. 

2.3 Nonlinear Response and radial distribution function 

In viscoelastic measurements on in vitro actin networks one observes strain 
hardening (Janmey et al., 1990), i.e. the system stiffens with increasing 



strain. This may either result from collective nonlinear effects or from the 
nonlinear response of the individual filaments. In the preceding section we 
have seen that the force coefficient obtained in linear response analysis for 
longitudinal deformation diverges in the limit of vanishing thermal fluctua- 
tions indicating that the regime of validity for linear response shrinks with 
increasing stiffness. Since the nonlinear response of a single filament may 
be obtained from the radial distribution function by integration, we discuss 
the latter first. 
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A central quantity for characterizing the conformations of polymers is the 
radial distribution function G(R; L) of the end-to-end vector R. For a freely 
jointed phantom chain (flexible polymer) it is known exactly and for many 
purposes well approximated by a simple Gaussian distribution. While rather 
flexible polymers can be described by corrections to the Gaussian behavior 
(Daniels, 1952), the distribution function of polymers which are shorter or 



comparable to their persistence length shows very different behavior. It is 
in good approximation given by 

G(R;L)«-A_/(|(i_^/L)), 

|exp[-7r^x] for X > 0.2 



where ^{x)=\ l/x — 2 



1 

4x 



(4) 

for X < 0.2 ^ ' 



and A/" is a normalization factor close to 1 ( |Wilhelm and Frey, 199(| ). This 
result is valid for L ^ 2^p, x ^ 0.5 and d = 2> where d is the dimension 
of space. A similar expression exists for d = 2. As can be seen in Fig. ^, 
the maximum weight of the distribution shifts towards full stretching as the 
stiffness of the chain is increased to finally approach the 5-distribution like 
shape required for the rigid rod. 

The radial distribution function is a quantity directly accessible to ex- 
periment since fluorescence microscopy has made it possible to observe the 
configurations of thermally fluctuating biopolymers ( |Gittes et al., 1993 ; Kas 



et al., 1993; Ott et al., 1993| ). Comparing the observed distribution func- 



tions with the theoretical prediction for d = 2 is both a test of the validity 
of the wormlike chain model for actual biopolymers as well as a sensitive 
method to determine the persistence length which is the only fit parameter. 
It should be noted here that the determination of persistence length e.g. 



of actin is still an actively discussed subject ( pupuis et al., 1996| ; [Wiggin! 
et al., 1997] ). 



A very interesting possibility would be to attach two or more markers 
(e.g., small fluorescent beads) permanently to single strands of polymers and 
to observe the distribution function of the marker separation. This would 
eliminate all the experimental difficulties associated with the determination 
of the polymer contour. Note that in contrast to existing methods of analysis 
it is not necessary to know the length of polymer between two markers; it 
can be extracted from the observed distribution functions along with ip by 
introducing L as a second fit parameter. 

The nonlinear response of the polymer to extending or compressing forces 
can be obtained from the radial distribution function by integration. The 
result (Fig. |2|) is in agreement with and provides the transition between the 
previously known limits of linear response and very strong extending forces 
(e.g., ( Marko and Siggia, 1995D ). For compressional forces, a pronounced de- 
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Figure 2: Left: End-to-end distribution function of a semiflexible polymer (nu- 
merical results). Note that with increasing stiffness of the polymer there is a pro- 
nounced crossover from a Gaussian to a completely non-Gaussian from with the 
weight of the distribution shifting towards full stretching. Right: The mean end-to- 
end distance i? as a function of a force applied between the ends (f = — /R/|R|). 
The step at positive (i.e. compressive) forces can be viewed as a remnant of the 
Euler instability. 

crease of differential stiffness around the classical critical force fc = ktt^/L^ 
can be understood as a remnant of the Euler instability. For filaments 
slightly shorter than their persistence length the influence of this instability 
region extends up to and beyond the point of zero force corresponding to 
the maximum in the linear response coefficient for £p ~ L (see Fig. [l|). For 
large compressions beyond the instability, the force-extension-relation cal- 
culated from the distribution function is only in qualitative agreement with 
numerical results because of the restricted validity of Eq. ^for x ^ 1. 

2.4 Dynamic structure factor 

The short time dynamics of polymers is most effectively measured by scat- 
tering techniques. The dynamic structure factor of a semiflexible polymers 
can be derived analytically in the limit of strong length scale separation 
a <^ X <^ ip, L ('weakly bending rod' approximation) and A < S^m- Here 
we have introduced the symbols a, A and for the lateral diameter of 
the polymer, the scattering wavelength and the mesh size of the network, 
respectively. The length scale separation guarantees that the decay of the 
structure factor is due to the internal dynamics of single filaments. For 
A ~ the structure factor decays in two steps due to the fast internal and 
to the slower collective modes, respectively. For the analysis we will restrict 
ourselves to rather dilute sol utions and refer the r eader for a more complete 
treatment to the literature ( |Kroy and Frey, 1997 ). 

The most important result of the calculations is that the time decay 
of the dynamic structure factor shows two regimes: the simple exponential 
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initial decay regime for t <^ t = C.^/nk'^ with a decay rate 

and the stretched exponential decay 

5(k,t) cx exp l^-^M(^^t)!^ , = kBTkl/ty^Q^ (6) 

for t ^ T. The effective transversal friction coefficient per length, ac- 
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Figure 3: Results from dynamic light scattering experiments. Left: Correction to 
the classical prediction 7^°^ ~ for the initial decay rate of the dynamic structure 
factor. The theoretical predictions for dilute solutions (solid line) and semidilute 
solutions (dashed line) are compared with experimental data ( pchmidt et al., 1989 ) 
Also included is the prediction for Gaussian chains from Ref. ( poi and Edwards 
1986). Right: Fit of theoretical dynamic structure factor to experimental data for 



k = 24.2 /im-i (Gotter et al., 1996) 



counts for the hydrodynamic damping of the undulations on length scales 
comparable to the scattering wavelength. For polymers with ~ A mainly 
the first regime, Eq. ||, is observed. The microscopic lateral diameter a of 
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the polymer enters the calculation as a microscopic cutoff and can be de- 
termined from a measurement of the initial decay of the dynamic structure 
factor. The method is easy to apply and has been shown to provide rea- 
sonably accurate results (see Fig. For polymers with ip X mainly the 
stretched exponential decay regime, Eq. ^ is observed. The decay rate is 
determined by the persistence length of the molecule, which is also read- 
ily extracted from experimental data. Some applications of the methods to 
dynamic light scattering with actin are shown in Fig. ^. 

Our analytical results for the initial decay rate and the dynamic expo- 
nent for semiflexible polymers suggest that known deviations of the dynamic 
exponent for more flexible polymers from its classical value z = 3 are most 
likely due to the local semiflexible structure of these molecules, and that the 
above analysis is therefore of some relevance also to scattering from flexible 
polymers. The reason is that the singular (oc 1/r) hydrodynamic interaction 
favors short distances. The possibility to determine the microscopic lateral 
diameter a of the polymer from the initial decay becomes important when 
the effects of bundling and side-binding are investigated, which can give rise 
to an effective thickening of the molecules. 



3 Many-chain properties 

As we have mentioned in the introduction, we are still far from understand- 
ing on a microscopic basis the macroscopic viscoelastic properties of solu- 
tions and gels of semiflexible polymers. In the following we will review some 
idealized concepts we have developed to model certain aspects of this compli- 
cated behavior and compare our predictions with some recent experiments. 



3.1 Entanglement transition 

Upon increasing the polymer concentration c at fixed polymer length there 
is a critical concentration c* above which a plateau develops in the frequency 
dependence of the storage modulus and the modulus increases steeply. The 
same phenomenon is observed at a critical length when the polymer length 
is varied at fixed concentration. It is called the entanglement transition. An 
experimental observation of the transition for short actin filaments (L ~ 1.5 
^m) as a function of concentration is shown in Fig. ^ (left). Adapting 



concepts which have been developed for flexible polymers (Kavassalis and 
Noolandi, 1989| ) we have proposed a mean field des cription for the entan 



glement transition in semiflexible polymer solutions ( Kroy and Frey, 1996|) . 



The basic idea is that the ends of a polymer are less efficient in confining 
other polymers than internal parts of the polymer. As a consequence the 
critical concentration c* is related to the overlap concentration c by some 
universal number C which may be interpreted as an effective coordination 



Statistical mechanics of semiHexible polymers 



11 



G [Pa] 



0.45 
0.35 
0.25 
0.15 
0.05 



12 3 4 

c [mg/ml] 



10 



G [Pa] 



10- 



□ o 



0.1 



□ ^ 
□ o 



c [mg/ml] 



Figure 4: Left: The entanglement transition for sliort rod-lilce actin filaments 
(L w 1.5 /im) of varying concentration. Right: Concentration dependence of the 
plateau modulus for pure actin (□) and actin with a small amount of gelsolin (o) 
(iHinner et al., 1997|). 



number and is a measure of the mean number of neighboring polymers nec- 
essary to confine the lateral motion of the test polymer. The prediction for 
an ideal gas of rods, oc kBT/^'^gL, is also shown in Fig. ^ From the 
fit we have determined the coordination number C, which compares very 
well with the result for flexible polymer solutions. Away from the transi- 
tion the effective mesh size reduces to the ordinary mesh size ^eflf oc £,rn- 
With oc c~^/^ (Schmidt et al., 1989) this implies oc c, which is in- 
deed observed for a regime of concentrations within the semidilute phase. 
However, at higher concentrations or lengths the slope of G^{c) increases, 
as can be seen from the bottom plot in Fig. ^, which shows the measured 
plateau modulus G^ for pure actin and actin with a small amount of gelsolin 
('^AG = 6000 : 1 corresponding to an average actin filament length of 16 fim). 
The description of the polymers as rods, which neglects internal modes of 
the individual polymers, is no longer appropriate. NOTE: Since the time 
of writing this manuscript, new experimental and theoretical investigations 
([Hinner et al., 1997 ) have suggested another, probably superior, interpreta- 
tion of the entanglement transition, which does not involve the concept of 
the coordination number. 



3.2 Plateau modulus 

To explain the observed elastic modulus on the basis of the elastic properties 
of the individual chains and the architecture of the network is in general a 
quite difficult task. Some of the theoretical work is based on a tube picture 
as depicted schematically in Fig. ^. The effect of the network surrounding 
an arbitrary test polymer is represented by a cylindrical cage of diameter 
d. The main response of the polymer to various deformations of this tube 
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Figure 5: Some of the theoretical work is based on a tube picture as depicted 
schematically here. The effect of the network surrounding an arbitrary test polymer 
is represented by a cylindrical cage of diameter d. Different scaling laws for d have 
been proposed. 



is assumed to arise from distortions of undulations of wavelength Lg, with 
Lg oc cP£p ( |Odijk, 1983 ). It is far from obvious, whether the macroscopic 
elasticity can be explained solely in term s of the compressibility of the tubes 
( pdijk, 1983 ; Isambert and Maggs, 1996 ), or whether contributions from fil- 
ament bending ([Kroy and Frey, 1996 ; Satcher and Dewey, 1996) or buckling 
( pVIacKintosh et al., 1995| ) are also important. The theoretical approaches 
differ quite significantly in their assumptions on how forces are transmitted 
through the network and which type of deformation of a single polymer gives 
the dominant contribution to the network elasticity. Due to those differences 
in the basic assumptions the predictions for the shear modulus of the net- 
work also differ. Our present data are reasonably consistent with entropy 
arguments assuming oc ksT/^'^Le, which leads to oc c'^/^ for d oc 

and G^ oc c'^/^ for d oc £m^/^y^ ( [Isambert and Maggs, 1996 ), respectively. 



Above 0.4 mg/ml our data also agree with a scaling G'^ oc c^^^, which has 
been derived by a scaling argument based on Fig. ^ for the case that the 
macroscopic elastic response is mainly due to filament bending ( |Kroy and 



Frey, 19961) . It should also be noted that these concentrations are close to 
the critical volume fraction 3a /L (a ~ 10 nm is the lateral diameter of an 
actin filament) for the nematic transition of rigid rods. 



3.3 Terminal relaxation 

At frequencies below the plateau regime the elastic response decreases 
and the polymer solution starts to flow. The corresponding time scale is 
called the terminal relaxation time. It can be determined from the mea- 
sured plateau modulus and the zero shear rate viscosity using the relation 
rjQ ~ G^Tr or directly from G'{uj). 

Intuitively, the mechanism for the terminal relaxation is obvious from the 
tube picture (see Fig. P) described above: viscous relaxation only occurs, 
when the polymers have time to leave their tube-like cages by reptation, 
i.e., by Brownian motion along their axis. The reptation model, which was 
originally formulated for flexible polymers, was adapted to the semiflexible 
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Figure 6: Left: Illustration of the generalized reptation picture for semiflexi- 
ble polymer solutions. Right: Terminal relaxation time above the entanglement 



transition (Hinncr ct al., 1997) 



case ( Odijk, 1983 ). Odijk ( pdijk, 1983| ) estimated the disengagement time 



Td for a semiflexible chain diffusing out of its tube. However, the data for 
Tr presented in Fig. ^ (bottom) are not in accord with his result for r^. 
The dependence of the observed terminal relaxation time on polymer 
length L is substantially weaker than predicted for r^, even in the stiff limit 
where = £pL/AD\\ oc (dot-dashed line in Fig. |6|), = ksT /2TTriL 
being the longitudinal diffusion coefficient of the chain in the free draining 
approximation. Instead, the solid line in Fig. ^ corresponds to the scaling 

law Tr OC 

A tentative interpretation of the data can be given in terms of a semiflex- 
ible polymer diffusing along a strictly one dimensional path; i.e., not being 
allowed to choose between infinitely many new directions at its ends. This 
situation is schematically depicted in the upper part of Fig. |6|. The charac- 
teristic decay t ime for self-correlat ions of the end-to-end vector (R(t)R) is 
then given by ( Hinner et al., 1997| ) 



m = L%/D\\ {Rr ^{L + 2i,Y/AD\\ . (7) 

This presents an upper bound for the terminal relaxation time within the 
tube model. As seen from the dashed line in Fig. ^, tr (for Ip = 17 /im 
( pittes et al., 1993 ; Ott et al., 1993] )) is in fact by a factor of ten too large 



compared to the data but describes fairly well the length dependence of 
Tr- The restriction to one path implies a very slow decay of conformational 
correlations. An unusually slow decay of stress (the frequency dependence 
of G'{uj) is still less than linear in the measured frequency range) is indeed 
observed, but this might also in part be due to the broad length distribution 
of actin ( |Janmey et al., 1986| ). Clearly, further investigations are necessary 
to come to a better understanding of the terminal regime. 
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3.4 Stochastic network models 

As motivated in the introduction, we use a two dimensional toy model of 
crosslinked semiflexible polymer networks to study some of the fundamental 



problems in semiflexible network theory ( Wilhelm and Frey, 1998 ). Here 
we address the question whether longitudinal or transversal deformation 
dominate the network response and find an interesting crossover between the 
two extreme cases. Since such questions are already nontrivial for networks 
at T = (i.e. no fluctuations), we consider this most simple case first: The 
linear elastic response of a network of classical compressible and bendable 
rods. A more sophisticated treatment will include contour fluctuations of 
the filaments (but no steric interactions) by replacing the rods with elements 
having the appropriate linear or nonlinear force extension relations described 
above. 

Rods of unit length are randomly placed on a square piece of the plane. 
Wherever two sticks intersect, they are connected by a crosslink. Periodic 
boundary conditions are applied to the left and right sides of the square. 
A shearing horizontal displacement is enforced on all rods intersecting the 
upper or lower boundary of the square at the intersection points. The (lin- 
earized) elastic equations of the system are solved numerically for displace- 
ments and forces at all crosslink points. The shear modulus is obtained from 
the force needed to impose the displacements of the upper and lower bound- 
ary. Crosslinks are not allowed to stretch but do not fix the angle between 
the intersecting rods. The results reported here change quantitatively but 
not qualitatively when crosslinks do fix the intersection angles. 

The parameters of this system are the density p of rods per unit area, the 
compressional and the bending stiffness of the rods. We choose to express 
the latter by the linear response force constants kcomp and fcbend for a rod 
of unit length. Geometrically, the relative size of the two force constants is 
controlled by the aspect ratio a = a/L oi the rods where a is the rod radius. 

While the system can show no elastic response (in the limit of infinite 
size) for densities below the percolation threshold ~ 5.72 ( pike and Sea- 



^er, 1974), the shear modulus G grows very rapidly with increasing p for 
P ^ Pc- For sufficiently slender rods G displays unusual power law behavior 
in p even for denisties far above the percolation threshold. 

To address the question whether the elasticity of a random stiff poly- 
mer network is dominated by transverse or by longitudinal deformations of 
the filaments, one can study the dependence of the shear modulus on the 
two force constants. We keep A;bend fixed and increase A:comp from values 
corresponding to a = 0.15 (short, thick rod) to values corresponding to 
a = 1 X 10~^ (long slender rod) for different system densities (see Fig. ^). 
We observe that beyond a certain point the shear modulus ceases to depend 
on /ccomp; indicating that the elasticity is dominated by bending modes for 
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Figure 7: Shear modulus of a crosslinked two dimensional random network of 
rods of unit length for different densities p and different aspect ratios a plotted 
over the relative size fccomp/^bend of compressional and bending force constants. 
Units were chosen in such a way that the bending stiffness k and hence fcbend were 

constant (k = 1). 

thin rods. Since the relevant scale for the system elasticity at higher den- 
sities is set not by the rod length but by the mesh size, the point of onset 
for this behavior shifts upwards with density. The dominance of bending 
modes in this region is confirmed by the observation that almost all of the 
energy stored in the deformed network is accounted for by the transverse 
deformation of the rods. Even in the region were the shear modulus does 
depend on A^comp, a substantial part of the elastic energy is stored in bending 
deformations. 
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